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We show that intensity of a wave created by a source embedded inside 
a three-dimensional disordered medium exhibits a non-universal space-time 
correlation which depends explicitly on the short-distance properties of disor- 
der, source size, and dynamics of disorder in the immediate neighborhood of 
the source. This correlation has an infinite spatial range and is long-ranged 
in time. We suggest that a technique of "diffuse microscopy" might be devel- 
oped employing spatially-selective sensitivity of the considered correlation to 
the disorder properties. 



Tremendous progress has been made during the last years in the understanding of mul- 
tiple scattering of waves in disordered media in so-called mesoscopic regime.0 This regime 
is commonly defined by 



where A is the wavelength, i is the mean free path, L is the sample size (or source-detector 
separation in the case of infinite medium), and Labs is the absorption length. In particular, 
considerable attention has been paid to the study of the correlation function of intensity 
fluctuationsEFi 
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I. INTRODUCTION 
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where /(r, t) = \ilj{r,t)f is the intensity of the wave field '?/'(r, t), SI{r,t) = I{r,t) — (J(r, t)), 
angular brackets (■ ■ ■) denote ensemble averaging, and the random fields ip{r,t), I{r,t) are 
assumed to be stationary. In the present paper, we are mainly interested in the correlation 
of intensities in two distant points separated by a distance superior to the mean free path: 
Ar > i. 

Several contributions to C can be identified.!'! Short-range (in both space and time) 
correlations arise from interference between the waves arriving at the points r and r + 
Ar. Weak long-range correlations are due to the interaction of diffusing modes in the 
medium. Finally, tiny infinite-range correlations manifest themselves in the phenomenon of 
universal conductance fiuctuations,! and have been recently observed for optical waves, ixhe 
listed contributions to the correlation function are commonly denoted by Ci, C2 and C3, 
respectively. It is important to emphasize that the magnitudes of these correlations scale 
with some dimensionless parameter and are not sensitive to the short-distance properties of 
disorder provided that |r — ro| , |r + Ar — ro| > For the case of plane wave transmission 
through a slab, such a parameter is the dimensionless conductance g ~ Nl/L, where N ~ 
k^A is the number of transport channels, fco = 27r/A is the wavenumber, and A is the area 
of the slab surface. For 1, which is a typical experimental situation, one findsi Ci ~ 1, 
C2 ~ g^^, C3 ~ g^^. For a medium of infinite extent, the relevant parameter is the product 
k^t Ci ~ 1, C2 ~ (A;o^)~^, C3 ~ (A;o^)~^5 where k^l ^ 1 is assumed. In all cases, the only 
medium-specific quantity determining the magnitude of the correlation is the mean free path 
£, and thus we call the C7i, C2, and C3 correlations universal as they are independent of the 
particular features of disorder, and only depend on the global parameter ^. 

Recently, a new contribution to the intensity correlation in random media (Co) has been 
conjectured.! This contribution originates from the scattering of wave near its source. In the 
case of a point source embedded in an infinite three-dimensional medium with fiuctuations 
of the dielectric constant obeying white- noise Gaussian statistics, for Ar > ^ and r = 0, it 
has been foundi 
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C«(r.Ar.O) = ^. (3) 

The physical process giving rise to the Co correlation consists in the scattering of the wave in 
the immediate neighborhood of its source (more precisely, at a distance of order A from the 
source), which produces a secondary wave propagating to the distant points r and r + Ar by 
diffusion. It turns out that the secondary waves created in such a way are correlated, and 
hence the described physical process establishes correlation between intensities at different 
points of the medium. It is important to emphasize that this correlation has an infinite 
spatial range. 

A fundamental difference between Cq and other types of correlation is that Cq should be 
sensitive to the short-distance details of disorder (e.g., to its correlation length).! Moreover, 
as this contribution to the correlation function is due to the scattering in the vicinity of the 
source, one could expect that Cq should depend on the geometrical parameters of the source 
and, in particular, on its size. Finally, for moving scatterers and r 7^ 0, Cq should also 
be sensitive to the dynamics of disordered medium in the immediate neighborhood of the 
source. It follows then that the Cq contribution to the correlation function is non-universal. 

In the present paper, we show that the Co correlation function does exhibit a strongly 
non-universal behavior. First, we study the dependence of Co on the correlation length of 
disorder i^. It is shown that the Co correlation function depends explicitly on koie, and, 
more precisely, that it decreases linearly with koie for kois ^ 1- The sensitivity of the Co 
correlation function to the correlation length of disorder might be used for sensitive specific 
measurements of Next, we study the effect of the source size a on the Co correlation. We 
show that Co decreases as a is increased. The decrease of Co with a is linear for koa <^ 1 and 
becomes power-law [Co oc (koo)^^] for koo ^ 1. Finally, we compute the space-time intensity 
correlation function Co(r, Ar,r). We show that this correlation function is sensitive to the 
motion of scatterers in the immediate neighborhood of the source. In the case of disordered 
medium consisting of scatterers undergoing Brownian motion, Cq decays linearly with \Jt/to 
as r -C To and exhibits a power-law behavior [Cq oc {t/to)^^^"^] for r ^ tq, with tq being a 



3 



single-scattering correlation time. The decay of Co with time is exclusively determined by 
the motion of scatterers located at a distance of order A from the source, and is not sensitive 
to the dynamics of the rest of the medium. Measurements of Co could therefore provide a 
sensitive spatially-selective probe of scatterer motion in disordered media. 



II. SPATIAL CORRELATION 

Following Ref. ^, we consider an infinite three-dimensional disordered medium with a 
monochromatic source of scalar waves located at tq. We assume that the source is spherically 
symmetric and has a radius a. To treat the problem of non-zero a, one needs to specify the 
boundary conditions for the wave field at the surface |r — ro| = a. We assume that waves are 
transmitted through the region of space occupied by the source without being scattered.lll 
Then for |r — ro| > a the wave field '4'{r,t) obeys the scalar wave equation with a point 
source at Tq: 

{ + kl [1 + 6e{r, t) + tv]}^{r, t) = 6{r - ro), (4) 

where 6e{r,t) is the fluctuating part of the dielectric constant [6e{r,t) = for |r — ro| < a], 
and ?7 is a positive infinitesimal ensuring the uniqueness of solution of Eq. @. The wave 
equation can be written in the form of Eq. (^) only if 6e{r, t) does not change significantly 
at the time scale of (cfco)"^? with c being the speed of propagation in the average medium: 
d/dt[6e{r,t)] < cko. 

Let us first consider the purely spatial correlation function C(r, Ar, 0). The time evo- 
lution of 6e{r,t) is of no importance in this case, and the statistical properties of fe(r, t) 
are completely specified by its average value and covariance. Without loss of generality, we 
consider 66{r, t) which obeys Gaussian statistics with zero average and correlation length i^: 

(fe(r,t))=0, (5) 



S,(Ar) = (fe(r,t)fe(r + Ar,t)) = cr^xp ( ) , (6) 
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where |r — ro| , |r + Ar — ro| > a. In the present paper we restrict ourselves to weak scat- 
tering {klllal ^ 1) and short correlation length (^^ <^ A). The first assumption allows us 
to benefit from a standard diagrammatic techniqueJi in the first order Born approximation, 
while the second one permits to not distinguish between the scattering and the transport 
mean free paths. The case of ^ A calls for a more sophisticated treatment which is beyond 
the scope of this paper. For k'^i'^a'^ ^ 1 and is ^ A, the mean free path is i = 47r/(-u/cQ) 
with u = TT^^'^a'^il. In the following we assume that any variation of implies corresponding 
change of so that i remains unchanged. 

As shown in Ref. in a weak scattering limit {koi ^ 1), the leading contribution to 
{SI{r,t)SI{r + Ar,t + t)) can be represented by a diagram depicted in Fig. |I], where the 
ladder propagators are L(ri,r2) = |ri — r2|), and is a vertex to be calculated with 
account for the short-distance properties of disorder and the size of the source. This vertex 
is a short-range object and can be replaced by a number. We show the leading contributions 
to V in Fig. 0. Note that the diagrams (c) and (d) contribute to the vertex only if 7^ or 
a 7^ and therefore they have not been counted in Ref. where only the case of = a = 
has been considered. 

The calculation of the diagrams of Fig. ^ is performed similarly to that of Ref. |^. For 
the diagram (b) of Fig. |^ we obtain, for instance. 



V^'^ = k', (J-^ I dh, I dh,G{r, - ro)G\u - ro)/*(ri - ro)/(r4 - ro)i?,(ri, r4). 



(7) 



where integrations cover the whole space except the interior of a sphere of radius a centered 
at ro. 



G(r - ro) = -—r- r exp 

47r |r — ro| 

is the average Green function of Eq. (11), and 



ikn ) |r — rnl 



(8) 



sin(/co |r - rol) / |r-ro|\ , . 

t,|r-ro| ""H"^^]- 

The derivation of the two above expressions is performed under important assumptions of 
koi ^ and a <^ i. The second assumption allows us to neglect the absence of the dielectric 



constant fluctuations inside the region of space occupied by the source, which simphfles 
greatly the derivations. 

Resuhs similar to Eq. (|^) are also obtained for the three other diagrams of Fig. |^. The 
diagram (c), for example, is given by the same integral as in the Eq. (|^ with the integrand 
G*(ri — ro)G'*(r4 — ro)/(ri — ro)/(r4 — ro)-B£(ri, r4). Next, diagrams (a) and (d) are obtained 
from (b) and (c) by a complex conjugation. Finally, taking the sum of all the four diagrams, 
performing angular integrals, and introducing Xi = k^ri and X2 = k^r^ as new variables of 
integration, we get 

V = + . . . + \/W = ^^F{koie, koa), (10) 



where 



X < exp 



1 I j_ f°° j_ sin(2a;i) sin(2x2) 

rf. \2'\\ 

(11) 
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F{a,(3) = I dxi / dx2 



T[-^l'' a Jf3 Jp Xi X2 
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^Xi - X2^ ^ 



a 



exp 



Xi+ X2\' 



a 
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The rest of the calculation is straightforward. The diagram of Fig. |l] yields 
{61{r,t)dl{r + Ar, t + r)) = \/L(ro, r)L(ro, r + Ar)[£/(47r)]2, where the factor [i/iAn)]^ ac- 
counts for the two vertices connecting the diffusion ladders to the points r and r + Ar. As 
the average intensity is (/(r, t)) = [£/ (47r)]^L(ro, r),i the contribution of the diagram of Fig. 
ID to the correlation function deflned by Eq. (H) becomes 

Co(r,Ar,0) = -^F(fco4,M- (12) 

koi 

The result is recovered for = a = 0, as F(0, 0) = 1. Hence, the factor F(A;o4, koo) 
describes the effects of source size a and disorder correlation length on the Co correlation. 
In general, F(fco4, ^0^^) is a rather complicated function of its arguments. Below we consider 
separately the following two cases: («) a ^ <^ A and (^^) 4 <^ a, X. 

(i). If a <^ <^ A, one can set the lower hmits of integrations in Eq. (|1^) equal to 
zero. Then Eq. (|TT]) can be simplifled, although the simpliflcation requires some care, as 
the integrand of Eq. (|ll]) does not have flnite derivatives at a = and cannot be therefore 



expanded in a power series near this point. Instead, we represent the expression in the curved 
brackets of Eq. ( pAj ) as 2exp[— (x^ + xl^lo?\ sinh(2a;iX2/Q;^) and then expand sinh(- ■ ■) in 
Taylor series. Carrying out integrations, we obtain^ 

0) = -a f + 1; 3/2; -a^), (13) 

TT ^or(n + 3/2) 

where ■ ■) is the confluent hypergeometric function. While the series in Eq. (|TB|) con- 
verges sufficiently fast for a > 1, the convergence slows down as a becomes less than unity. 



As we consider <^ A, we are only interested in -F(a, 0) for a ^ 1. Although Eq. (|T3D 
can be evaluated numerically even for a < 1, this equation is of little help if one wants to 
get some qualitative insight into the small-a behavior of F[a^ 0). To get such an insight, we 
introduce t\^2 = xi^2 /oi as new variables of integration in Eq. (|TT]). This yields 

2 1 /-oo roo sin(2ati) sin(2at2) 



2 i ' 
F{a,0) = —j^- dtj dt2- 
TT'^/^ a Jo Jo 



Tc-^'^ a Jo Jo ti t2 

X exp (-tl - tj) sinh(2tit2). (14) 



Assuming small a, we replace sin(2at2) by 2Q;t2 and then perform the integrals over t2 and 
ti. Finally, we keep in a final result only the terms linear in a. This gives 

F(a,0)^l-^a, a<l. (15) 



A plot of F{a,0) obtained by a numerical evaluation of Eq. (|T3|) is presented in Fig. y 
by a solid line. The decay of F is linear for a < 1, as predicted by Eq. (|l^). The part of the 
plot corresponding to a > 1 is meaningless for a moment, as we have assumed ^ A. This 
part of the plot will be discussed and interpreted in the next section, where the space-time 
correlation function is considered. 

As follows from Eq. (|I5|) and Fig. |^, the Cq correlation decreases with increase of the 
disorder correlation length i^, and the characteristic scale of this decrease is set by A. This 
is physically clear as Cq originates from scattering at a distance of order A from the source, 
which is only possible if the wave "identifies" its environment as disordered already after a 
distance of order A. Such an identification implies that 6e{r,t) changes significantly at the 



scale of A, or, in other words, that the correlation length of the fluctuations of Se{r, t) is 
shorter than A. On the basis of this reasoning, we expect that Cq should vanish for > X, 
while for ~ A the behavior of Cq may be sensitive to a particular model of disorder. 
In contrast, the linear initial decay (|T5D of the Co correlation function with the disorder 
correlation length seems to be rather general and independent of the model of disorder. 
We have verified that the linear law (^) holds (with an apparent change of numerical factors) 
not only for the Gaussian covariance of Se{r,t) [Eq. (|^)], but also for the exponential one 
-Be(Ar) = (T^ exp(— Ar/4), as well as for -Be(Ar) = a'^Il{Ar / i,.) , where Il{x) = 1 for x < 1 
and n(x) = for a; > 1. 

An essential result of our analysis is that the Cq correlation function appears to be 
sensitive to the correlation length of the dielectric function fluctuations, which is a short- 
distance parameter of disorder. We remind that all the previously studied types of correlation 
(Ci, C2, C3) are not sensitive to such details of the scattering medium, and only depend on 
the mean free path i. We conclude then that the sensitivity to the short-distance properties 
of disordered medium is a distinctive feature of the Cq correlation function. Measurements 
of Cq could thus provide unambiguous information on the disorder correlation length. It is 
rather surprising that such information is not lost as a result of multiple scattering, as one 
would expect, but determines the long-range spatial correlation in the multiple-scattering 
speckle pattern. 

Another important outcome of the presented consideration is that the correlation length 



is entering into Eq. (|T2D is that characteristic for the disorder in the immediate neighborhood 
of the source (more precisely, at a distance of order A from it). The information about the 
rest of the medium is only contained in the mean free path i. Measurements of the Cq 
correlation could therefore be used for "diffuse microscopy" of disordered media with spatial 
resolution of order A. A possibility of such a microscopy stems from the enhanced sensitivity 
of the multiple-scattering speckle pattern to the properties of disorder in the vicinity of the 
source. 

(ii). In the case of £e -C a, A, we can set 4 = and consider the limit a — > of 



Eq. ([TTD . The exponentials in the curved brackets of Eq. (|Tl|) become 7i^/'^aS(xi — X2) and 
n^^'^a5{xi + X2), respectively, which yields 



2 

F(0, koa) = 1 

TT 



Si(4koa) - 



sin (2koa) 
2koa 



1 — (4/7r)A;oa, koa -C 1, 

(16) 

l/{27ikoa), koa 3> 1, 
where Si(- ■ ■) is the sine integral, and the second line of the equation shows simplified versions 
of the first one in two limiting cases. Eq. ([T6|) is plotted in Fig. |^ by a dashed line. As 
k^a <C 1, the decay of F(0, fcoa) is linear, similarly to that of F(fco4)0) (see the solid line 
in Fig. ^). When koa becomes larger than unity, F{0,koa) decreases as {koa)~^. This 
indicates that although the small size of the source is important for the magnitude of the Cq 
correlation, the correlation process leading to Cq survives even for extended (a 7^ 0) sources. 

The above analysis holds for Ar > i. If Ar < i, an additional contribution to the Co 
correlation function arises due to scattering in the neighborhood of the observation points. In 
a complete analogy with the term given by Eq. ([l^), this second contribution is determined 
by the short-distance properties of disorder in the vicinity of the observation points r and 
r + Ar. For Ar = 0, the contributions owing to scattering near the source and near the 
observation point become equal and their sum determines a correction to the second moment 
of intensity fluctuations (5/(r,t)^). Hence, the most of the conclusions drawn concerning 
the Co correlation are also true (with obvious insignificant modifications) for ((5/(r, t)^). 

III. SPACE-TIME CORRELATION 

We now turn to the case of the space-time correlation function, allowing r 7^ in Eq. 
(0). It is obvious from general considerations that the correlation function of intensity of 
a multiple-scattered wave should decrease with r once the disordered medium is subject 
to some non-periodic dynamics. This decrease was first discovered for the Ci correlation 
function,011 and recently has been also studied for C2 (see Ref. |13D and C3 (see Ref. 
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As the role of the disorder correlation length (.^ has been studied in the previous section, 
and in order to emphasize the effect of the medium dynamics, we assume = in the present 
section. For the sake of concreteness, we consider the disordered medium as a suspension 
of point-like scatterers which diffuse with a diffusion constant Db- The relevant correlation 
function of the dielectric constant fluctuations is theni'0 

5.(Ar, r) = (fc(r, t)fe(r + Ar, t + r)) = j^^^^^, exp (^-^ j . (17) 

Calculation of the diagrams of Figs. |l|, |^ is performed very similar to the case of r = 0. 
We finally obtain 



TT 



Co(r, Ar, r) = —F [^Jt/t,, haj , (18) 

where F{- ■ ■) is given by Eq. (|1^), and tq = (4A;qDs)^^ is the single-scattering correlation 
time. To simplify the further analysis, we assume a -C A, {tDbY^"^ which allows us to 
extend the lower limits of integrations in Eq. (|ri|) to zero, and to reduce Eq. ( pTD to Eq. 

An approximate expression for F{a, 0) in the limit of small a has been found in the 
previous section [see Eq. (|T5|)]. To analyze the behavior of F(q;,0) for a 1, we take the 
large-a limit of Eqs. ([Tl]), ([13|): 



F{a,Q)-^—^, a>l. (19) 

Hence, the asymptotic behavior of the space-time correlation function can be summarized 
as follows: 



TT 



l-(4M')v/^, r«ro. 



Co(r,Ar,r)^ — X ' (20) 

l(l/7r3/2)(r/ro)-3/^ r » tq. 

These approximate results are well supported by the numerical evaluation of Eq. (|13|) which 



yields the solid line of Fig. ^ Linear initial decay of F (and, consequently, of Co) with \Jt/tq 
is followed by a power-law long-time tail, as predicted by Eq. ([191) . 
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It is instructive to compare the initial decay of the Co correlation function with that of 
Ci, which for the considered case of a point source in an infinite medium and for Ar = 
readsEi 

C,(..0..).l-2^/f, (21) 

For Ar 7^ 0, Ci decreases exponentially with Ar/£ and thus becomes negligible for Ar > £.i 
Although the functional forms of Co(r, Ar, r) and Ci(r, 0,r) for short r are the same, the 
typical correlation time tq corresponding to Co is much larger than the correlation time 
ro(£/ |r — ro|)^ corresponding to Ci, as |r — ro| ^ ^. This stems from the fact that the decay 
of Ci is due to multiple scattering events, each giving a contribution to the decorrelation 
process. The factor \r — tqI /£ in Eq. (^) reflects the increase of the effective number of 
scattering events with |r — ro|. In contrast, the decay of Cq with r originates from single 
scattering events in the neighborhood of the source. The following propagation of waves 
to the distant points r and r + Ar depicted by the ladder propagators in Fig. |I], does not 
lead to any additional decorrelation, and thus the characteristic decay time of Co(r, Ar, r) 
is insensitive to r — ro and r + Ar — rQ. 

It is worthwhile to note that there exists another important difference between the Cq 
and Ci correlation functions. Eq. (^) for Ci is derived under assumption of statistically 
homogeneous scatterer dynamics throughout the whole space [the characteristic time Tq in 
Eq. (^) is assumed to be the same for all parts of the medium]. If, for instance, the 
scatterer motion in some region of space differs from that in the rest of the medium, the 
difference is reflected directly on the decay of Ci. This property of Ci allows its use for 
imaging of scatterer dynamics in turbid me diali. Moreover, space-resolved measurements of 
Ci permit to distinguish between Brownian and directed (flow) motion of scatterers inside 
some bounded region of space, and to image scatterer flows in multiple-scattering media.0 
In contrast, as the decay of Cq is only sensitive to the motion of scatterers in the vicinity 
of the source, any modifications of scatterer dynamics in other parts of the medium do not 
affect the Cq correlation [the characteristic time tq in Eqs. (|18|) , (pOf) is that corresponding 
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to the immediate neighborhood of the source]. The Cq correlation function is not therefore 
suitable for imaging of the bulk of disordered media. As a compensation of this drawback, 
Co has, however, an important advantage as compared to Ci. The advantage is the spatially- 
selective sensitivity to the scatterer motion in the immediate neighborhood of the source. 
This sensitivity opens up a possibility of "dynamic diffuse microscopy" of scatterer motion 
in multiple-scattering media. 

The presented treatment of the space-time correlation is valid for Ar > i. If Ar < 
i, an additional contribution to the correlation function appears due to scattering in the 
neighborhood of the observation points, just like in the case of purely spatial correlation. 
For Ar = 0, the sum of the two contributions yields a correction to the Ci correlation 
function [Eq. ([21|)]. A possible existence of such a correction has been suggested by Scheffold 
et a/., El although they considered a different experimental geometry. This correction is 
uniquely determined by the motion of scatterers in the neighborhoods of the source and the 
observation point. We emphasize that Cq <^ Ci for Ar = 0, while for Ar > i the Ci term 
disappears and the leading contribution to the correlation function is given by Cq, since 
C2 ~ (koi)-' « Co. 

IV. CONCLUSION 

In the present paper, we consider a recently suggested! novel type of intensity correlation 
in disordered media (Co). We show that the Cq correlation function exhibits a strongly non- 
universal behavior in the sense that it is sensitive to the short-distance properties of disorder 
and geometrical parameters (e.g., size) of the localized source of waves. In particular, we find 
that Co depends explicitly on the correlation length of dielectric function fluctuations and 
the size of the source a. More precisely, the Co correlation is damped down as either £e or a is 
increased. As a function of koi^, Co exhibits linear decay for koie ^1- As a function of koa, 
the Co correlation function also decays linearly for koa <^ 1, and Co oc {koa)^^ for koa ^ 1. 
The sensitivity of the Co correlation to the short-distance properties of disorder (e.g., to i^) 
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is its distinctive feature, since all other types of correlation commonly considered for waves 
in disordered media (Ci, C2, C3) are insensitive to the details of scattering medium, and 
only depend on the mean free path, which is an integral parameter of disorder. 

Next, we calculate the space-time correlation function of intensity fluctuations 
Co(r, Ar, r). We show that this correlation function is only sensitive to the scatterer motion 
in the immediate neighborhood of the source, while being insensitive to the dynamics of 
the rest of disordered medium. Co(r, Ar, r) is independent of r, Ar and decays hnearly as 
a function of \Jt/tq for r <S tq (tq is the single-scattering correlation time). For large r 
(r ^ To), the decays of Co is power-law: Cq oc (t/tq)"^/^. 

Finally, we suggest that the spatially-selective sensitivity of the Cq correlation function 
to the properties and dynamics of disordered medium in the neighborhood of the source 
(and detector for Ar < t) may be used to develop a sort of "diffuse microscopy" of multiple- 
scattering media. 
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FIGURE CAPTIONS 

FIG. 1. The leading diagram contributing to the correlation function {51{r, t)SI{r + Ar, t + t)). 
L's denote the ladder propagators, the left propagator corresponds to time t, and the right one — 
to t + T. Solid lines with arrows denote retarded (with arrows directed from ro to r or r + Ar) 
and advanced (with arrows directed from r or r + Ar to tq) average Green functions, G and G* , 
respectively. V is the vertex located at tq. Diagrams contributing to V are shown in Fig. 2. 

FIG. 2. The diagrams contributing to the vertex V. Solid lines with arrows directed from tq to 
some other point denote retarded average Green functions G, solid lines with arrows directed to 
ro denote advanced average Green functions G* . In each diagram, the left pair of Green functions 
corresponds to time t, and the right pair — to t + r. Crosses indicate scattering (a factor is 
associated with each cross), dashed lines denote correlation function i?e(Ar, r) of the fluctuating 
part of the dielectric constant. Integrations are assumed over ri, r2, r^, and r4. ro is the source 
position. 

FIG. 3. Function F{a,P) describing the effects of non-zero correlation length i^, size of the 
source a, and time delay r on the Cq correlation function. The solid line is F{a, 0), and the dashed 
line is F [0,(5). The attenuation factor which multiplies the Cq value found for = a = and r = 
is obtained by setting a = kot^ (disorder with non-zero correlation length, koi^ <^1), a = (t/to)^^^ 
(scatterers undergoing Brownian motion), or /3 = koa (non-zero size of the source). 
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